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Abstract. The framework of anisotropic hydrodynamics is used in 3+1 dimensions to analyze be- 
havior of matter produced in ultra-relativistic heavy-ion collisions. The model predictions for the 
hadronic transverse-momentum spectra, directed and elliptic flows, and the HBT radii are presented. 
We show that the effects of the initial anisotropy of pressure may be compensated by appropriate ad- 
justment of the initial energy density. In this way, the final hadronic observables become insensitive 
to the early stage dynamics and the early thermalization puzzle may be circumvented. 
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INTRODUCTION 

Soft-hadronic observables measured in the ultra-relativistic heavy-ion experiments are 
well reproduced by perfect-fluid hydrodynamics [1] or by viscous hydrodynamics with a 
small viscosity to entropy ratio [2, 3, 4, 5, 6, 7], Nevertheless, the use of such approaches 
at the very early stages of the collisions encounters important conceptual difficulties. 
Thermalization times shorter than a fraction of a fermi (used quite often in the perfect- 
fluid approaches) cannot be explained within the microscopic models of the collisions. 
On the other hand, viscous hydrodynamics is based on an implicit assumption that one 
can make an expansion around the isotropic background. If the shear correction is large, 
a new framework is needed. This has triggered developments of the reorganization of 
viscous hydrodynamics in which one incorporates the possibility of large momentum- 
space anisotropics into the leading order of the approximation. Such a new framework 
has been introduced in [8, 9, 10, 11, 12, 13, 14] and we refer to it below as to the 
anisotropic hydrodynamics. 

ANISOTROPIC HYDRODYNAMICS 

Anisotropic hydrodynamics is based on the equations 

d^ v = 0, (1) 
= E, (2) 



which express the energy-momentum conservation and entropy production laws. The 
energy-momentum tensor T^ v has the structure 

T iiv = (g +jP± ) [/M[/V _ p±gM v _ ( P± _ P|| )yMyV 5 (3) 

where is the longitudinal pressure and P± is the transverse pressure. In the limit 
^|| = P± = P, Eq. (3) reproduces the energy-momentum tensor of the perfect fluid. 
Similarly, the entropy production law (2) is reduced to the entropy conservation law, 
if we take E = 0. The four- vector [/^ is the four-velocity of the fluid, while is the 
four- vector whose spatial part is parallel to the beam (z) axis. In the general case, we use 
the parameterizations JJ^ = (wocosh#, w x , u y ,uosinh'&) and = (sinh#, 0,0, cosh 
where u x and u y are the transverse components of the four-velocity field and ■& is the 
longitudinal fluid rapidity. The entropy flux is defined by the formula 

a^ = aU^, (4) 

where o is the non-equilibrium entropy density. 

One can show [8] that instead of and P± it is more convenient to use the entropy 
density o and the anisotropy parameter x as two independent variables (one may use 
the approximation P\\/P± = x~ 3 ^ 4 ). Similarly to standard hydrodynamics with vanish- 
ing baryon chemical potential, the energy density £, the entropy density a, and the 
anisotropy parameter x are related through the generalized equation of state £ = e(a, x) . 
In our model we use the following ansatz [10] 

s(x,a) = e qgp ((7)r(x), (5) 
P±(x,a) = P q gp(o) [r(x) + 3xr'(x) 
P\\(x,a) = P qgp (o)[r(x)-6xr'(x) 



where the functions £ qgp and P qgp define the realistic QCD equation of state constructed 
in Ref. [15]. The function r(x) characterizes properties of the fluid which exhibits the 
pressure anisotropy x. Here we use the form introduced in [8] 
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r(x) = — 



jcarctan-\A— 1 
1 + r== 



(6) 



In the isotropic case x = 1, r(l) = 1, /(l) = 0, and Eq. (5) is reduced to the standard 
equation of state used in [15]. 

The function E in Eq. (2) defines the entropy source. We use the form proposed in [8] 

na , x) = Hz£f!L, ( 7) 

where the time-scale parameter T eq controls the rate of equilibration 1 . In the limit of 
small anisotropy Eq. (7) is consistent with the quadratic form of the entropy production 



1 In this work we use the constant value T eq = 1 fm. In Refs. [9, 11, 13] the medium dependent T eq was 
used, which was inversely proportional to the typical transverse momentum scale in the system. If the 
constant value of T eq is used, the system approaches the perfect fluid behavior for T ^> T eq . 
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FIGURE 1. (Color online) The pseudorapidity distribution of charged particles for 
T e q = 1.0 fm. The results are obtained for xo = 100 (dashed blue lines) and xq = 0.032 (dotted 
green lines), and for the two centrality classes: c = — 6% (b = 2.48 fm) and c = 15 — 25% 
(b = 6.4 fm). The model results are compared to experimental data from PHOBOS [19] (red 
dots) and BRAHMS [20] (green squares). 



in the Israel-Stewart theory. Far from equilibrium, hints for the form of E are lacking, 
although we may expect some suggestions from the AdS/CFT correspondence [16]. 
Thus, for large anisotropies the formula (7) should be treated as an assumption defining 
the dynamics of the system. 

INITIAL CONDITIONS AND FREEZE-OUT 

In the general 3+1 case we have to solve Eqs. (1) and (2) for o, x, u x , u y , and ■&, which 
depend on t, xj_ = (jc, y) , and 7] (t is the proper time and rj is the space-time rapidity) . We 
fix the initial starting time to To = 0.25 fm. Similarly to other hydrodynamic calculations, 
we assume that there is no initial transverse flow, w x (to,x^,ti) = u y {lQ,x^,r\) = 
and that the initial longitudinal rapidity of the fluid is equal to space-time rapidity, 
$(To,x^,T7) = r/. In this text we present the results for two scenarios: i) the initial 
source is strongly oblate in the momentum space, jc(to,xj_, T}) = 100, and ii) the source 
is prolate in momentum space, *(to,xj_, 77) = 0.032; the latter value is chosen because 
r(100) = r(0.032). 

The initial entropy density profile has the form 

cr (T|,Xi.) = cr(To,T7,x 1 ) = e gq ^[£ip(&,T7,x 1 )], (8) 

where b is the impact parameter, and p(b, T],XjJ is the normalized density of sources, 
p(b,r],xj_) — p(b : 7],x±)/p (0,0,0), for details see [14]. The quantity £[ is the initial 
energy density at the center of the system created in the most central collisions. Its value 
is fixed by the measured multiplicity, separately for two different physical scenarios 
considered in this paper. We use £[ = 48.8 GeV/fm 3 and 80.1 GeV/fm 3 forjt = 100 and 
xq = 0.032, respectively. 
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FIGURE 2. (Color online) Transverse momentum spectra of 7T + 's for the centrality class 
c = — 5% (b = 2.26 fm) and for different rapidity windows; xq = 100 (dashed blue lines) 
and = 0.032 (dotted green lines), and T eq = 1.0 fm. The model results are compared to 
the experimental data from BRAHMS [21]. The spectra for different rapidity windows are 
successively rescaled down by factor 0.1. 



The evolution is determined by the hydrodynamic equations until the entropy density 
drops to Of = 1.79 fm~ 3 , which for x = 1 corresponds to the temperature Tf = 150 
MeV. According to the single-freeze-out scenario, at this moment the abundances and 
momenta of particles are expected to be fixed. The processes of particle production and 
decays of unstable resonances are described by using THERMINATOR [17, 18], which 
applies the Cooper-Frye formalism to generate hadrons on the freeze-out hypersurface. 



RESULTS 

Our results are shown and compared to the RHIC data (Au+Au collisions at the highest 
beam energy = 200 GeV) in Figs. 1-5. The model results describing the psudo- 

rapidity distributions, transverse-momentum spectra, the elliptic and directed flow co- 
efficients, and the HBT radii have been obtained with different initial anisotropies of 
pressure; xq = 100 (dashed blue lines) and xq = 0.032 (dotted green lines). In the all 
considered cases we find a good agreement between the model results and the data, for 
more results see [26] . Moreover, we find that the results obtained with different initial 
anisotropies are practically the same. This is so because we have adjusted the initial en- 
ergy density separately for two different values of xq. A larger (smaller) initial energy 
density is used for the initially prolate (oblate) system. 

Our results indicate that the final hadronic observables are not sensitive to the early 
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FIGURE 3. (Color online) Transverse-momentum dependence of the elliptic flow coefficient 
V2 of 7T + + K + calculated for c = 20 - 40% (b = 7.84 fm) at midrapidity and for T eq = 1.0 fm; 
xo = 100 (dashed blue lines) and xq = 0.032 (dotted green lines). The results are compared to 
the PHENIX Collaboration data (red dots) [22]. The presented errors are statistical. Horizontal 
bars denote the pj bins. 
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FIGURE 4. (Color online) Pseudorapidity dependence of the directed flow of charged particles 
for the centrality bin c = 5 — 40% and T eq = 1.0 fm; xq = 100 (dashed blue lines) and xq = 0.032 
(dotted green lines). The results are compared to the experimental data from STAR (red dots) 
[23] and PHOBOS (green squares) [24]. 



anisotropy of pressure. The flows are being built during the whole time evolution of the 
system, hence the relatively short early anisotropic stage does not influence the results. 
In our opinion, the insensitivity of the hadronic observables helps us to circumvent the 
early thermalization puzzle. 
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FIGURE 5. (Color online) The HBT correlation radius 7? out of positive pions as a function of 
the total transverse momentum of the pair for the centrality c = — 5%, T e q = 1-0 fm, and for 
two values of the initial anisotropy: xq = 100 (dashed blue lines) and xq = 0.032 (dotted green 
lines). The results are compared to experimental data from STAR Collaboration (red dots) [25]. 
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